Transport Bifurcation in a Rotating Tokamak Plasma 
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The effect of flow shear on turbulent transport in tokamaks is studied numerically in the experi- 
mentally relevant limit of zero magnetic shear. It is found that the plasma is linearly stable for all 
non-zero flow shear values, but that subcritical turbulence can be sustained nonlinearly at a wide 
range of temperature gradients. Flow shear increases the nonlinear temperature gradient threshold 
for turbulence but also increases the sensitivity of the heat flux to changes in the temperature gra- 
dient, except over a small range near the threshold where the sensitivity is decreased. A bifurcation 
in the equilibrium gradients is found: for a given input of heat, it is possible, by varying the applied 
torque, to trigger a transition to significantly higher temperature and flow gradients. 



Introduction. Turbulent transport of heat is a major 
obstacle to the development of a successful fusion de- 
vice. Turbulence powered by microinstabilitics such as 
the ion temperature gradient (ITG) instability rapidly 
transports heat out of the plasma, limiting the temper- 
ature gradient that can be sustained by a given input 
of heat, and thus the temperature that can be reached 
at the core of the plasma. The problem is exacerbated 
by the strong dependence of the turbulent amplitudes on 
the driving gradients, which in general keep the gradient 
not far above the critical threshold for the onset of the 
instability, a phenomenon known as stiff transport 

Experimental evidence 0, Q suggests that a sheared 
flow in the plasma greatly improves the situation; such 
a flow can significantly reduce turbulent fluxes for given 
values of the driving gradients. In some cases a large 
enough shear can help quench the turbulence altogether 
d, 0] . Flow shear may also reduce the sensitive depen- 
dence of heat flux upon temperature gradient (the " stiff- 
ness") 0. Several of these results have been confirmed 
in numerical simulations [7-ll|, but the picture remains 
incomplete. Furthermore, while in simulations one may 
specify the flow shear, in experiment the control param- 
eters are the input of heat and momentum. A set of 
simulations would ideally show not only that a large flow 
shear is beneficial, but also how it may be achieved. 

A recent paper has demonstrated the basic prop- 
erties of turbulence in the Cyclone Base Case regime 
(concentric circular flux surfaces with q = 1.4, t/Rq = 
0.18, RoVii/vthi — 0.01, where q is the magnetic safety 
factor, r the half diameter of the flux surface measured in 
the midplane, i?o the major radius at the magnetic axis, 
Vthi the ion thermal velocity and va the ion-ion collision 
frequency) at a finite value of magnetic shear. It was 
found that the ITG-driven turbulence was quenched at 
sufficiently high flow shear, but the ITG was replaced as 
the driver of the turbulence by the parallel velocity gradi- 
ent (PVG) [3, At large flow shears, the system was 
linearly stable, but strong subcritical PVG-driven tur- 
bulent transport could be sustained at sufficiently large 



temperature gradient. 

In this Letter we consider the case where the shear of 
the magnetic field is zero, which experimental observa- 
tions have indicated is favourable for quenching turbu- 
lence 0, 0, [l^ . It is found that the plasma is now lin- 
early stable for all non-zero values of flow shear, and that 
there is a much larger range of flow shear and ITG values 
where the turbulence is completely quenched. It is shown 
that because of this (and when neoclassical transport is 
taken into account) there is a steady-state bifurcation in 
the flow and temperature gradients, for certain values of 
input heat and applied torque [l3|. A positive feedback 
between the suppression of turbulence and the input of 
momentum can cause a one-way jump in both gradients. 

Theoretical Framework. The work reported in this 
letter follows on from that carried out in [l2| , which con- 
tains a more detailed exposition of the model used. All 
simulations have been carried out using the code GS2 
[l8| . which solves the local non-linear gyrokinetic equa- 
tion in the presence of a sheared toroidal flow Ruj [l9l.[2o|. 
where R is the major radius of the tokamak and lu is the 
toroidal angular velocity. The flow is ordered as smaller 
than the ion thermal velocity in a Mach-number expan- 
sion: Ruj ^ Mvthi where pi/R ^ M <^ 1 and pi is 
the ion Larmor radius. In order to ensure that the ef- 
fect of the flow shear is retained, the gradient of the 
flow is then ordered as the inverse of the Mach number 
(dlnw/dlnr ^ 1/-^), so that the flow shear is of or- 
der the fluctuation frequency and the particle streaming 
rate: 7b = {r / q)duj / dr vthi/R- Using this expansion, 
the effects of the sheared toroidal flow are included in 
GS2 [13, by adding a time dependence to the ra- 
dial wavenumbers and a drive term associated with the 
PVG. All simulations reported here are electrostatic with 
a modified Boltzmann electron response. Typical resolu- 
tion was 32 X 64 X 14 X 4 X 8 (poloidal, radial, parallel, 
pitch angle, energy). 

The gyrokinetic ordering is used _to_ close the moment 
equations in the transport model 
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FIG. 1: Linear behaviour: the time evolution of the heat flux, 
normalised to its initial value, for different values of 7b at 
R^i/Lt — 11. Growth is transient for all non-zero flow shears, 
and the length of the period of transient growth decreases with 
flow shear. Dashed curve is subcritically stable (no turbulence 
sustained nonlinearly) . Inset: Inverse transient growth time, 
and the number of exponentiations the heat flux undergoes 
before switching to decay, both vs. flow shear. 
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FIG. 2: Turbulent heat flux (a) and toroidal angular momen- 
tum flux (b) vs. flow shear for different values of Rq/Lt- For 
lower values of Rq/Lt, there is a range of flow shears where 
the turbulence is fully quenched. 



and lit, respectively, can then be calculated from the 
perturbed ion distribution function Sf as the flux sur- 
face averages of the convection of those quantities by the 
fluctuating E x B velocity [ioj . 

Subcritical Turbulence. With the flow shear equal 
to zero, the fluctuation amplitude grows exponentially 
with time (Fig. [ij. As the flow shear increases from zero, 
however, the growth becomes transient, and switches to 
decay after a time r-y which decreases with increasing 
(inset to Fig. [T]). This is in qualitative agreement with 
recent theoretical [l^, and numerical [l^l results at a fl- 
nite value of magnetic shear, except that at zero magnetic 
shear we observe that there are no growing eigenmodes 
for any non-zero value of 7^. 

With this in mind, it might be expected that there 
would be no turbulence for je > 0. However, linear 
instability is not necessary to sustain nonlinear turbu- 
lence in a rotating plasma; in fact, the transient growth 
caused by the ITG or PVG drive is sufficient to give 
rise to subcritical turbulence in simulations initialised 
with sufficient-amplitude noise (or fully developed tur- 
bulence). Thus, there is turbulence, but it is subcritical 
for all finite flow shears (in contrast with the finite mag- 
netic shear case [l^). 

Heat Flux. Considering the dependence of turbulent 
heat flux on the flow shear. Fig. Elja), we observe 
that for lower values of the temperature gradient, viz. 
Ro/Lt ^ 11.5, where Lt is the temperature gradient 
scale length, Qt decreases smoothly to zero with increas- 
ing flow shear. The turbulence is then fully quenched 
for a range of flow shears, but then the PVG drive be- 
comes so large that turbulence is reignited. For larger 
values, Rq/Lt > 11.5, the turbulence is never quenched; 
Qt merely decreases and rises again, which echoes results 



12l |. We stress, however, that turbulence is quenched 



at lower values of flow shear and higher values of Rq / Lt 
than in 1^, and so there is a much wider range of pa- 
rameter space where the flow shear is large enough to 
quench the ITG-drivcn turbulence, but not large enough 
to drive PVG turbulence. This more favourable regime is 
what enables the transport bifurcations described below. 

Examining the dependence of the heat flux on Rq/Lt, 
Fig. El shows that at je < 0.4, there are two nonlin- 
ear thresholds (Fig. [3l^c)), which both increase with flow 
shear. Below the first threshold turbulence is completely 
quenched; between the first and second thresholds Qt in- 
creases slowly, and above the second it rapidly rises to 
values similar to the case without flow shear. Thus, flow 
shear signiflcantly increases the overall temperature gra- 
dient required for turbulence, reduces the transport stiff- 
ness dramatically between the flrst and second thresholds 
(i.e. at low values of Qt), but increases the stiffness above 
the second threshold (Fig. |3{^d)). For je > 0.8 there is 
only one threshold, which decreases to with increas- 
ing flow shear, (Fig. [3ljb)), as the PVG drive increases. 
Stiffness is low but all the fluxes are very large. 

Momentum Flux. The story of the toroidal angular 
momentum flux lit is quite simple. Defining the turbu- 
lent viscosity as vt = {Jlt/lE){r/rLimiR^q), the turbu- 
lent heat diffusivity as xt = {Qt/ {Ro/ Lt)) (Ro/in-iTi)) 
(where rui, Ui and Ti are the mass, density and tem- 
perature of the ions), the turbulent Prandtl number, 
Pr ~ i^t/xt, was in the range 1.0-1.8 in all our simula- 
tions, with only a very weak dependence on the gradients. 
In other words Ilt/Qt cjc je/{Ro/ Lt)', the turbulence 
transports heat and momentum in equal proportions. 

Fig. mb) details this behaviour: lit rises with je, 
reaches a local maximum, drops, then rises again as the 
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FIG. 3: Turbulent heat flux vs. Ro/Lt for (a) -^e < 0.8 and 
(b) 7b > 0.8. (c) Nonlinear turbulence thresholds vs. flow 
shear. For each value oi Rq/Lt, the Ist-threshold curve shows 
the 2 values of 7_b at which the turbulence is quenched and 
then rekindled (cf. FigEJ. For Rq/Lt > 11.5 or 7b > 1.8, the 
turbulence is always present, (d) Proflle stiffness vs. the flow 
shear. Low Qt refers to fluxes between the flrst and second 
thresholds, High Qt to fluxes above both thresholds. 



PVG starts to add significantly to the ITG turbulence 
drive. As with the heat flux, there is a window of zero 
turbulent transport for Rq/Lt ^ 11.5. The extent of this 
window in and Rq/Lt is shown in Fig. [UJc). 

Quenching Turbulence. At this point we abandon the 
use of jE and Rq/Lt as independent variables, since the 
actual control parameters in a fusion device are the total 
rates of input of heat and momentum by external sources. 
If we assume a steady state, the rates of input of these 
quantities are equal to their outgoing fluxes, i.e. Q and 
n. In a practical manner, a question may be posed: in 
choosing Q and 11, how can the temperature gradient, 
and hence the temperature at the core of the plasma, be 
maximised? Since the input of infinite amounts of heat 
is not possible, the turbulence must be quenched. 

If turbulent transport is reduced, coUisional transport 
becomes important. For tokamaks, a quantative the- 
ory of this transport, known as neoclassical theory, ex- 
ists [2^. For the case of circular concentric fiux sur- 
faces in the banana regime (va <^ qRo/vthi) studied 
here, the neoclassical thermal diffusivity and viscosity 
are Xn ^ OM{R/r)^^^'^\^ pfvi^ and Q.lq^pfvu- 
So the total heat fiux Q ~ Qt + Qn where Q„ = 
Xn{Ro/ LT)niTi/ Rq, and the total toroidal angular mo- 
mentum fiux n = IIj + Tln where n„ = v^'yEniniiR^q/r. 
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FIG. 4: (a) The ratio of the total momentum flux 11 to the 
total heat flux Q vs. flow shear for a constant value of Q. An 
increase in applied torque (i.e., an increase in II/Q) at point 
A will cause a transition to point B. (b) Total heat flux Q vs. 
Rq/Lt for different constant values of II/Q. The points A 
and B on both graphs correspond to the same states. Since 
neither Q nor 11 can be specifled for a simulation (with the 
exception of 11 = 0), the contours of constant Q and II/Q were 
interpolated from a large number of data points using radial 
basis functions with a linear kernel [Sj. Also plotted are the 
neoclassical contributions to II/Q and Q. The contours in 
(b) do not intercept the neoclassical line but curl round and 
asymptote to it, tending back to the origin, and thus point B 
is in fact on the same contour as point A (see text and [251]'). 
This feature cannot be shown as the contours are too closely 
spaced for interpolation near the neoclassical line. 



The essential point is that the neoclassical Prandtl 
number, Vn/Xn — 0.01, is much smaller than the tur- 
bulent Prandtl number: turbulence is much more effec- 
tive at transporting momentum than collisions alone. It 
should be noted that provided this is satisfied, the qual- 
itative results of the next section do not depend on the 
exact values of the neoclassical transport coefficients. 

Triggering a Transition. Let us now consider what 
happens when we increase II/Q at constant Q. Fig. IH^a) 
shows a curve where the heat flux is held constant at 
Q = 2.6 UiTiVthiPi / Rq- As n/Q is increased, initially, 
because Q is being held constant, the suppression of the 
turbulence caused by an increase in 7^; causes a corre- 
sponding increase in Rq/Lt, which restores the turbu- 
lence to its former levels. However, at the point where 
Rq/Lt becomes so large that Qt ^ Qn, this negative 
feedback is broken, because Qn does not depend on ^e- 
The turbulence is greatly reduced and yet the heat flux 
remains unchanged because the bulk of the heat is be- 
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ing transported neoclassically. The same is not true of 
the momentum, owing to the much lower neoclassical 
Prandtl number: because of the reduction in turbulent 
amplitudes, the transport of momentum drops dramati- 
cally. This causes a large increase in the flow shear gra- 
dient, which rekindles turbulence which then transports 
momentum once again. In other words, when the mag- 
nitude of n/Q is increased above 0.06 Ro/vthi, there is a 
transition from point A to point B on Fig. |4l^a). 

Fig. lUJb) shows the same transition on the {Q, Rq/Lt) 
plane where contours of constant H/Q are plotted. 7^ in- 
creases along these contours from high Q and low Rq / Lt 
to low Q and high Rq/Lt- As increases, initially 
Q drops rapidly because of turbulence suppression by 
flow shear, until Qn starts to become significant. At this 
point, because H/Q ~ IIt/{Qt + Qn), must increase 
again to keep H/Q constant. As flow shear increases 
further, however, Ht/Qt becomes so large (Fig. [2]) that 
the turbulence amplitude, and hence Qt, must start to 
decrease again to maintain a constant H/Q. Thus the 
curve asymptotes to the neoclassical line, and the sys- 
tem ends up in a state where the heat transport is nearly 
all neoclassical, but the momentum transport is nearly 
all turbulent. Points A and B correspond to points A 
and B onHlJa). Point B is not on the neoclassical line 
in Fig. m^b), but corresponds to a (slightly) turbulent 
state (as can be seen from Fig. |l{a)), where the large 
7e means that lit is significant even though Qt <^ Qn- 
The location of point B on|4l^b) was calculated from|4l^a) 
using the relation xt ^ vt- During the transition Rq/Lt 
jumps from 7.4 to 10.4. Including the original suppres- 
sion, flow shear has enabled a total jump (at constant Q) 
from Rq/Lt ^ 4.5 at U/Q = (point C on Fig. Hb)), 
to Rq/Lt ^ 10.4 at U/Q ~ 0.06. 

Conclusions. In summary, we have shown that al- 
though the plasma is linearly stable for all finite values 
of flow shear, transient growth is sufficient to allow sub- 
critical turbulence to be sustained nonlinearly. At low 
values, flow shear reduces transport in two ways: by dra- 
matically increasing the threshold temperature gradient 
required to drive turbulence, and, over a small range of 
temperature gradients, by reducing the strength of the 
dependence of the fluxes on the temperature gradient 
(the stiffness). At the top of that small range there is 
a second threshold, above which the fluxes rapidly rise 
to levels seen with no flow shear. High values of flow 
shear can in fact increase the transport. 

Perhaps more importantly, we have discovered a tran- 
sition to a higher-gradient regime. The transition oc- 
curs when either the input of heat has been reduced or 
the input of momentum increased to the point where the 
bulk of the heat is transported neoclassically and the ion 
temperature gradient is no longer driving sufficient tur- 
bulence to transport the toroidal angular momentum. A 
positive feedback loop starts where the build up of the 
velocity gradient reduces the turbulence further until it 



becomes sufficient to drive turbulent transport via the 
PVG instability. A new stable equilibrium is reached 
with much higher temperature and flow gradients, and 
where the heat transport is nearly neoclassical. 

The speciflc numbers associated with this transition 
arc likely to depend significantly on the plasma configu- 
ration; for example, the value of q may have a large effect 
0. Extending this work to other regimes, in particular 
more experimentally realistic configurations, would thus 
be of interest. However, qualitatively we have shown not 
only that there exist equilibrium states with equal fiuxes 
yet different gradients, but that, with the right condi- 
tions, less favourable regimes can automatically transi- 
tion to more favourable ones with higher gradients. 
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